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1. Introduction 
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The wide-angle scattering of X-rays by anisotropic powders of nano-layered crys- 
tallites (nano-stacks) is addressed. Assuming that the orientation distribution prob- 
ability function / of the nano-stacks only depends on the deviation of the crystal- 
lites' orientation from a fixed reference direction, we derive a relation providing 
/ from the dependence of a given diffraction peak's amplitude on the azimuthal 
angle. The method is applied to two systems of Na-fluorohectorite (NaFH) clay 
particles, using synchrotron radiation and a WAXS setup with a two-dimensional 
detector. In the first system, which consists of dry-pressed NaFH samples, the 
orientation distribution probability function corresponds to a classical uniaxial 
nematic order. The second system is observed in bundles of polarized NaFH par- 
ticles in silicon oil; in this case, the nanostacks have their directors on average in 
a plane normal to the reference direction, and / is a function of the angle between 
a nano-stack's director and that plane. In both cases, a suitable Maier-Saupe func- 
tion is obtained for the distributions, and the reference direction is determined 
with respect to the laboratory frame. The method only requires one scattering 
image. Besides, consistency can be checked by determining the orientation distri- 
bution from several diffraction peaks independently. 



The recording of a powder diffraction signal using a two- 
dimensional detector is a potentially simple and efficient tech- 
nique for determining characteristic structural length scales in 
crystals. The basis of the technique lies in that a scattering vol- 
ume containing an isotropic powder consisting of a large num- 
ber of particles with identical crystalline structure offers to the 
incoming beam all possible incidence orientations with respect 
to the crystallites. Consequently, all the diffraction peaks char- 
acteristic of the crystallite structure are visible in any of the 
azimuthal planes, and all azimuthal planes are equivalent: a sim- 
ple scan of the deviation angle, 2 9, in a given azimuthal plane, 
provides a spectrum containing peaks at all structural length 
scales characteristic of the crystallites, apart from those extinct 
due to form factor effects. The scattering image collected by a 
two-dimensional detector is isotropic, and its center of symme- 
try denotes the direction of the incident beam. Radial cuts of 
the images passing through this center provide identical one- 
dimensional spectra that can be averaged over all azimuthal 
angles in order to significantly increase the signal to noise ratio. 

It is well known that reconstruction of the three-dimensional 
crystal structure from powder data is made difficult by the over- 
lapping of reflections with similar diffraction angles. The inter- 
pretation is fairly easy in the case where only a dominant fea- 
ture of the structure is investigated. A particular case is that of 
nano-layered systems, where the powder is made of particles, 
or clusters of particles, that consist of a stack of identical quasi- 



2D crystallites. The Bragg reflections resulting from the set of 
Bragg planes associated to those stacks are the most intense and 
can be interpreted in quite a straightforward manner. A num- 
ber of natural (vermiculite, montmorillonite) and synthetic (flu- 
orohectorite) smectite clays have the properties of presenting, 
in the weakly-hydrated states, a microstructure where particles 
are platelet-shaped and nano-layered. The repeating crystallite 
is a lnm-thick silicate-platelet, and the stacks are held together 
by cations shared between adjacent platelets. The stacks can 
swell by intercalation of ions and molecules (water, in particu- 
lar), and powder diffraction has been used extensively to study 
that swelling, both on natural (Wada et al, 1990) and synthetic 
(da Silva et al, 2002) smectite clays. 

A wide range of natural or man-made materials are aggre- 
gates of crystallites with the same crystallographic structure. 
The distribution of the crystallite orientations inside the mate- 
rial are rarely isotropic, hence their X-ray diffraction signal is 
that of an anisotropic powder. Many examples of such poly- 
crystalline aggregates can be found among metals. The analy- 
sis from X-ray diffraction data of the preferred crystallographic 
orientations in those metals, i. e., their texture, provides infor- 
mation on the history of their deformation (Wenk, 1985), as it 
does for rocks, another type of polycrystalline materials with a 
larger degree of complexity (several types of crystals coexist). 
This is done by decomposing the distribution of the reflections' 
intensities as a function of the angular direction into a super- 
position of distributions resulting from single reflections {pole 
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figures), and subsequently inferring the orientation distribution 
probability (ODP) for the crystallites. Several general meth- 
ods (Wenk, 1985; Bunge, 1993) have been developed to tackle 
that potentially difficult (and sometimes, non-univoque) inver- 
sion. From the experimental point of view, 2D detectors have 
proved useful in such texture analyses, as they allow gathering 
of the data necessary to the pole figure analysis in a short time 
and without changing the sample position (Ischia et ah, 2005). 
See (Wenk & Grigull, 2003) about the practical applicability 
of those methods. Note that texture studies have recently found 
an interesting application in methods allowing to determine the 
complete 3D structure of a crystal from anisotropic powder 
diffraction data collected from samples with different textures 
(Lasocha & Schenk, 1997; Wessels et al, 1999). 

The study of preferential orientations in polycrystalline mate- 
rials consisting of nano-layered particles (quasi- ID crystals) 
such as the clays presented above has made little use of the gen- 
eral inversion technique for texture data, since there is basically 
no overlap of pole figures for those materials. In the case where 
the ODP possesses an axial symmetry, it is classically deter- 
mined experimentally by measuring a rocking-curve, in which 
the sample is rotated around an incidence angle 9 corresponding 
to a well known (OOfc) reflection, and the scattered intensity is 
measured at the deviation angle 2 9 as a function of 9 (Giintert 
& Cvikevich, 1964; Taylor & Norrish, 1966). In this paper, 
we present an alternative method for determining the orienta- 
tion distribution from the scattering data in the case in which 
it only depends on one characteristic angle. When the assem- 
bly of nano-stacks present a partial priviledged orientation, all 
azimuthal planes are not equivalent any more. Our method is 
based on monitoring the dependence of the intensity of a given 
reflection peak on the azimuthal angle. It can thus be carried 
out on a single two-dimensional wide angle scattering image. 
In other words, we monitor the anisotopry of that image at a 
given reflection. Note that monitoring of the intensity along 
Debye rings has previously been done in a similar manner for 
cubic and hexagonal materials (Puig-Molina et al, 2003; Wenk 
& Grigull, 2003). We do not address the influence of the ODP 
on the diffraction profiles in a given azimuthal plane, a problem 
to which much attention has been given in the past, for systems 
like pyrolytic graphite (Giintert & Cvikevich, 1964) and smec- 
tite clays (de Courville et al, 1979; Plancon, 1980). 

The paper is organized as follows. In section 2, we present the 
theoretical background of the method, for both geometries. The 
method is applied to data from two experiments corresponding 
to two geometrical configurations observed in systems of the 
synthetic smectite clay Na-fluorohectorite (in sections 3 and 4, 
respectively). In section 5, we discuss the interest of the method. 
Section 6 is the conclusion. 

2. Theoretical background 

We consider a clay sample placed at the origin O of the labo- 
ratory frame (0,x,y,z) (see Fig. 1). An incident X-ray beam 
parallel to x, and of wavelength A, is scattered by the sample, 
and the scattered intensity recorded in a detector plane yO'z, O' 
being the point where the direct beam would hit the detector in 



the absence of a beam stop behind the sample. The sample con- 
sists of a population of nano-layered crystallites consisting of 
~ 100 stacked identical layers. The orientation distribution / of 
the crystallite population with respect to the laboratory frame is 
assumed known and anisotropic, resulting in an anisotropic dis- 
tribution of the scattered intensities in the detector plane. The 
running point M on the detector is denoted by the direction of 
the line Ou passing through the points O and M; the correspond- 
ing angular frame consists of the scattering deviation angle 2 9 
and the azimuthal angle <f). 

2.1. Diffraction by an anisotropic distribution of nano-stacks 

The orientation of each crystallite's director, n, with respect 
to the fixed frame, is denoted by the angular frame (@,4>) (see 
Fig. 1), with © in the range [0; ir/2] and <t> in the range [— w; tt]. 
& is the angle between the crystallite's director and the plane 
yOz, while 4> is the azimuthal angular coordinate of the direc- 
tor. 

Let us first consider a single crystallite whose stacking planes 
make with the incident beam an angle 9, corresponding to the 
Bragg condition for the stack. This implies that it mostly scat- 
ters along the direction Ou, i.e, that the azimuthal direction for 
scattering is 4> = 3>, and the deviation angle is 2 9 = 2 © (see 
figure 2(a)). As the azimuth cf> describes the [— tt; tt] range, the 
trajectory on the unit sphere of the director for the crystallites 
that scatter along that azimuth describes a circle in the plane 
x = sin 9 (see Fig. 2(b)). This director is therefore defined in 
the frame (x,y,z) by the coordinates 

/ sine \ 
»2fl(0) -cos6» sin0 . (1) 
y cos 9 cos0 J 

Next, we assume that the angular distribution / only depends 
on an single angle a, which is the angle between a crystallite's 
director n and a reference orientation The distribution is 
normalized with respect to the total solid angle available to the 
description, i. e. 

/ / f(a) sin a da df3 = 1 , (2) 
Jo Jo 

where {a, (3) is a standard spherical angular frame. We assume 
that the size- and orientation- distribution of the particles are 
uncorrected, and that the beam defines a scattering volume that 
is large enough for the size distribution of the scatterers to be 
entirely contained in it. Hence, we assume that the amplitude of 
a diffraction peak in the direction (2 6,(f>) is proportional to the 
number of scatterers that satisfy the Bragg condition for that 
direction, and that the proportionality constant is independent 
of <f>. This also assumes that the scattering peak is well sepa- 
rated from other peaks along the 2 9 line. We denote he{<t>) the 
amplitude of the diffraction peak at deviation angle 2 9 and at a 
given azimuthal angle value <fi, and a2e(<fi) the angle between 
the reference director no and that of the crystallites that diffract 
in the (2 6, <f>) direction. The intensity he{(t>) A(f> along a length 
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Acj) of the angular scattering cone of half-opening 29 is propor- 
tional to the quantity f(a2e{4 > )) Aa(</>, A<t>), where the angular 
length Aa is that travelled by n on the unit sphere as <f> describes 
the range A<fi. The trajectory in Fig. 2(b) is travelled at constant 
speed when varying <p uniformly, as appears from the vector 
dnldcf) having a norm cos# (computed from Eq. (1)), and thus 
independent of (f>. Therefore, Aa only depends on A(f>, which 
yields: 

/ 2e (0)A</>cx/(a 2e (0))Aa(A0) . (3) 

Consequently, a uniform sampling of the azimuthal profile of 
a WAXS peak in the scattering figure is directly related to the 
orientation distribution for the scatterers, /. In particular, a con- 
stant profile corresponds to a uniform angular distribution of the 
platelets, i.e, the scattering figure of an isotropic clay powder is 
isotropic, as expected. 

The angle a 2 g(0) can be computed from its cosine which 
equals no • «2 (</>)• The orientation with respect to the labora- 
tory frame of the mean director iiq is defined by the angles ©o 
and 4>o, and thus its coordinates in the (x,y,z) frame are 



sin ©o 
«o — cos ©o sm *o 
cos©o cos4>o, 



Thus, 



otieicf) = acos [sin 9 sin ©o + cos 9 cos ©o cos((/> — <E>o]) . 

(5) 

For a uniform sampling of the azimuthal scattering intensity 

profile he{4>)> one can directly relate the function he{4>) and 
f(a) bringing together Eq. (3) and (5) according to 

he{4>) oc /(acos[sin6* sin ©o + cos 9 cos ©o cos(0 — <E>o)]) ■ 

(6) 

From a given functional form for the distribution /, with a 
given set of shape parameters (peak position, width, reference 
level), it is possible to fit a function in the form (6) to the exper- 
imental data, and infer both the distribution's shape parameters 
and the reference orientation «o(®o, •t'o)- The obtained param- 
eters are independent of the particular reflection (2 9 or q) cho- 
sen, which provides a interesting way of checking the results' 
consistency. 

2.2. Two particular geometries 

2.2.1. Uniaxial nematic configuration: In the classical uni- 
axial nematic geometry, the reference orientation is the mean 
orientation for the population of nano-stacks. The smaller a, 
the more crystallites have their orientation along n(a). The ori- 
entation distribution / is peaked in and monotically decreas- 
ing to its minimum f(n/2). An order parameter (de Gennes & 
Prost, 1993), defined as 



cos a 



l) f , 



(7) 



i.e S = 7T / (3 cos 2 a — 1 ) f(a) sin a da , (8) 
Jo 



quantifies the overall alignment of the crystallites. For S = 0, 
the population is completely isotropic, while for S — 1 all crys- 
tallites are perfectly aligned. The alignment can also be quan- 
tified in terms of the root mean square (RMS) width Wf of the 
distribution, according to 



(9) 



y = (a )f = 2tt a f(a) sin a da 



2.2.2. Uniaxial "anti-nematic" configuration: In the config- 
uration that we denote "anti-nematic", the crystallites have their 
directions on average in a plane perpendicular to the reference 
direction. In other words, the stacking planes of the crystallites 
contain the reference direction, on average. The distribution / 
is peaked around n/2, monotically increasing with a from its 
minimum value /(0). By analogy to the nematic description, an 
order parameter can be defined as well, according to: 

S = (3 cos 2 (| - - 2) = (3 sin 2 a - 2) f , (10) 

i.e S = 2-k I (3 sin 2 a — 2) f(a) sin a da, (11) 
Jo 



(4) and an RMS angular width according to: 



i.e w 



(!-)' 

,2 =2 7T / ^— — aj f (a) sin a da . 



(12) 
(13) 



2.3. Modelled azimuthal profiles 

Fig. 3 shows expected azimuthal profiles, computed from 
Eq. (6), in the case of a nematic geometry where the distribu- 
tion / is of the classical Maier-Saupe form observed in uni- 
axial nematic liquid crystals (Maier & Saupe, 1958; Maier & 
Saupe, 1959): 



f{a) oc exp(m cos a) 



(14) 



For a reflection corresponding to a particular value of the devi- 
ation angle 2 9 — 3.29°, two values of the distribution's angu- 
lar width Wf (15 and 30°), and three different orientations of 
the mean director iiq (©0 = —5,0 and 5°) have been used. 
The profiles are normalized to an integral value 1 over the 
whole azimuthal range. A profile corresponding to ©o = 5° 
for 29 — 6.54 is also shown. 

Fig. 3 shows clearly how the parameters describing the orien- 
tational ordering of the crystallites translate into the azimuthal 
dependence of a diffraction peak's intensity. The angular width 
of the orientation distribution controls the width of the profile's 
peaks. The azimuthal coordinate of the mean director, 4> () , only 
has the effect of translating the profile along the <f) scale, while 
@o is responsible for the curve's asymmetry with respect to the 
azimuth <p = <t>o + 90°, or, equivalently, the asymmetry of the 
2D diffractogram with respect to the plane xOy (see Fig. 1). In 
addition, the magnitude of this asymmetry depends on the value 
of 2 9: if there is asymmetry (i.e, if ©0 ^ 0), it shall be more 
pronounced for the higher order peaks of a given reflection. 

3. Application to a nematic geometry 
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3.1. Fluorohectorite, a synthetic smectite clay 

Fluorohectorite belongs to the family of smectite (or swelling 
2: 1) clays, whose base structural unit is lnm-thickphyllosilicate 
platelet consisting of two inverted silicate tetrahedral sheets that 
share their apical oxygens with a tetrahedral sheet sandwitched 
in between (Velde, 1992). Fluorohectorite has a chemical for- 
mula per unit cell X_ c -(Mg3_ JC Ll c )Si40ioF2. Substitutions of Li + 
for Mg 2+ in part of the fully occupied octahedral sheet sites 
are responsible for a negative surface charge along the platelets. 
This results here in a large surface charge, 1.2 e~(Kaviratna 
et al, 1996) (as opposed to 0.4 e~ for the synthetic smec- 
tite clay laponite), allowing the platelets to stack by sharing an 
intercalated cation X, for example Na + , Ni 2+ , Fe 3+ . The result- 
ing nano-stacks remain stable in suspension, even at low ionic 
strength of the solvent, and are therefore present both in water 
solutions, gels (DiMasi et al, 2001) and weakly-hydrated sam- 
ples (Knudsen et al, 2004). 

Wide-angle X-ray scattering of gel (DiMasi et al, 2001) 
and weakly-hydrated (da Silva et al, 2002) samples of Na- 
fluorohectorite (for which X=Na + ), have shown that the char- 
acteristic platelet separation is well defined, and varies step- 
wise as a function of the surrounding temperature and humidity 
(da Silva et al, 2003). Samples are anisotropic, but the width 
of the corresponding WAXS peaks allows inferring the typi- 
cal nano-stack thickness, which is of about 80 — 100 platelets. 
This thickness is significantly larger than thicknesses (< 20 
platelets) for which asymmetry and irrationality in the (OOfc) 
reflections are known to occur (Drits & Tchoubar, 1990), which 
makes the analysis of (OOfe) peaks in powder diffraction spectra 
of Na-fluorohectorite relatively easy. 

3.2. Dry-pressed samples of Na-fluorohectorite 

Dry-pressed samples of Na-fluorohectorite were prepared 
from raw fluorohectorite powder purchased from Corning Inc. 
(New York). The raw powder was carefully cation-exchanged 
so as to force the replacement of all intercalated cations by Na + 
(L0voll et al, 2005). The resulting wet gel was then suspended 
in deionized water again, and after sufficient stirring time, water 
was expelled by applying a uniaxial load in a hydrostatic press. 
The resulting dry assembly was then cut in strips of 4mm by 
47mm. 

The resulting samples are anisotropic, with a marked aver- 
age alignment of the clay crystallites parallel to the horizontal 
plane (Knudsen et al, 2004), as shown in the two-dimensional 
sketch of the sample geometry, in the insert to Fig. 4. A nematic 
description such as that presented in section 2.2.1, with a mean 
director n close to the direction of the applied uniaxial load, 
describes their geometry well. 

3.3. Experimental setup 

Two-dimensional diffractograms of dry-pressed samples 
were obtained on beamline BM01A at the European Syn- 
chrotron Radiation Facility (ESRF) in Grenoble (France). We 
used a focused beam with energy 17.46 keV corresponding to a 
wavelength 0.7lA. The scattering images were recorded using 
a two-dimensional MAR345 detector, with a resolution of 2300 



pixels across the diameter of the 345mm wide circular detec- 
tor. The distance between the sample and the detector was cal- 
ibrated to 375.25mm using a silicon powder standard sample 
(NBS640b). 

Fig. 4 shows a simplified sketch of the experimental geom- 
etry. The samples were lying with the average nano-stack 
director close to vertical. A sample cell allowing humidity- 
and temperature- control in the surrounding atmosphere is not 
shown. A more complete description of the sample environment 
can be found in (Meheust et al, 2005a). 

3.4. Image analysis 

The center part of a two-dimensional scattering image 
obtained from a dry-pressed Na-fluorohectorite sample at a tem- 
perature of 98°C is shown in Fig. 5(a). The image anisotropy 
appears clearly, with a much stronger intensity of some of 
the diffraction rings along the vertical direction. Some grainy- 
looking isotropic rings can also be seen, they are due to impu- 
rities (mostly mica). In what follows, we infer the orientation 
distribution for the clay crystallites from three rings, separately. 
These rings correspond to the three first orders of the same 
reflection. 

3.4.1. Method: The images were analyzed in the following 
way. After finding the image center (incident beam direction), 
we compute diffraction lines at azimuthal angle values <f> m , vary- 
ing <j> m between and 359° by steps of 1°. For a given value 
of (j> m , the diffraction line is obtained by integrating over (j> 
the image subpart corresponding to values of <f) in the range 
[4>m — S(p/2; <j) m + 5<p/2], where 5<p = 5°. Six of these powder 
diffraction profiles are shown in Fig. 6. They exhibit 6 peaks 
characteristic of the nano-layered structure of the crystallites: 
the 3 first orders for the population of crystallites with no water 
layer intercalated (OWL), and the orders 1, 2, and 4 for the pop- 
ulation of crystallites with one water layer intercalated. At this 
temperature, the 1WL state is unstable (da Silva et al, 2003); 
it is present in these pictures because the transition between the 
two hydration states is not finished yet. The analysis is carried 
out on the peaks 0WL(001), 0WL(002) and 0WL(003), which 
are the most intense, and are weakly perturbed by other peaks. 
The corresponding characteristic length scales are listed in table 
1 ; they are in good agreement with previous measurements (da 
Silva etal., 2002). 

For each of the three peaks 0WL(001), 0WL(002) and 
0WL(003), the peak's intensity was plotted as a function of (j>, 
and the obtained profiles were fitted with Eq. (6). The three nor- 
malized profiles are shown in Fig. 7; they have been translated 
vertically for clarity. As a functional form for the orientation 
distribution, we chose the Maier-Saupe function (Eq.( 14)). We 
also allowed for a constant background Iq to be accounted for. 
This resulted in fitting to the profiles a function in the form 
(using Eq. (6) and (14)) 



Iq + C exp 



m (sin 9 sin ©o + cos 6 cos ©o cos(</> — 4>o)) z 



(15) 

with the fitting parameters Iq, C, m, ©o and 4>q. Note that 
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Eq. (15) depends on the diffraction peak chosen, as the value 
for the half-deviation angle 9 is no parameter for the fit but is 
different for each of the three profiles. As each of the fit param- 
eters controls a particular shape aspect of the profiles (see sec- 
tion 2.3), we can derive estimates for the parameters Iq, C, m, 
©o and <t>o, using a few approximations (see appendix A), and 
use these estimates as initial parameters for the fits. The fitting 
procedure is fast and the fitted parameters do not depend on the 
initial guesses. 

Using the fitted model, we then derived the quantities char- 
acteristic of the orientational ordering in the samples, namely 
the order parameter S (Eq. (8)) and the angular RMS width w/ 
(Eq. (9)). 

3.4.2. Results: The fits are shown in Fig. 7 as solid lines, on 
top of the data points. As <j> moves away from the peak values 
<t>o and <t>o + t/2, and as the intensity goes down to the noise 
level of the diffraction profiles shown in Fig. 6, the significance 
of the intensity /(</>) in terms of orientation distribution of the 
profiles is somewhat diluted. Hence, the points in the fiat sec- 
tions outside the peaks in Fig. 7 are little relevant when fitting 
the model to the data. We therefore estimate an uncertainty on 
the fitted parameters from the fluctuations in the fit parameters 
obtained by varying the range of data points (between 30 and 
90 degrees from the peaks) considered for the fit. This uncer- 
tainty is considerably larger than the statistical error on the fitted 
parameters themselves, for a given fit. 

For each of the profiles, the fitted parameters and the derived 
quantities are shown in table 2. The results obtained from the 
different orders are consistent with each other within the uncer- 
tainties, except for the parameter ©o where a 0.5°discrepancy 
exists. This is however small, and the first order profile was 
expected to behave slightly differently from the others as the 
corresponding peak is not totally separated from the 1 WL(001) 
peaks (see Fig. 6). This profile also gives a value slightly 
different for the angular width, but consistent with the other 
orders within the error bars. We can thus estimate the aver- 
age orientation of the nano-stack population to (©o,<I>o) = 
(5.8°±0.5°; 176.3°±0.1°). This value results from a possible 
non-parallelism of the average crystallite orientation from the 
direction of the constraint applied during sample preparation, 
and from a possible misaligment of the sample (a few degrees 
at most) from the horizontal direction. 

4. Application to an "anti-nematic" configuration 

We now apply the technique to data obtained in wide angle X- 
ray scattering experiments from bundles of Na-fluorohectorite 
particle in silicon oil. 

4.1. Electorheological chains of ftuorohectorite particles 

Fluorohectorite samples were prepared by suspending 1.5% 
by weight of clay particles in silicon oil (Rothiterm Ml 50, with 
a viscosity of lOOcSt at 25°C). The hydration state of the parti- 
cles prior to suspension in the oil was 1WL. When an electric 
field above 500 V.mm -1 was applied, the suspended crystal- 
lites polarize, and the subsequent induced dipoles interact with 



each other, leading to aggregation of the clay crystallites in 
columnar bundles (Fossum et ai, 2005). This can be seen in 
the microscopy picture shown in Fig. 8. These bundles are on 
average parallel to the applied electric field. 

4.2. Experimental setup 

The orientation of the polarized crystallites inside the bundles 
was investigated using synchrotron wide angle X-ray scattering. 
A top view sketch of the experimental setup is shown in Fig. 9: 
the suspensions were contained in a sample cell in between two 
vertical copper electrodes separated by a 2mm gap. A 1 kV elec- 
tric field was applied between the electrodes. Horizontal bun- 
dles of polarized particles were observed to form within 20s 
after the field was switched on, and two-dimensional diffrac- 
tograms were recorded after formation of the bundles. The 
experiments were performed at the Swiss-Norwegian Beam- 
lines at ESRF, using the same beam energy and detector as pre- 
sented in section 3.3. 

4.3. Image analysis 

Fig. 5(b) displays the center of a scattering image obtained 
after the chains have settled in the electric field. This WAXS 
data contains an significant background, with a broad ring due 
to diffuse scattering from the oil, but the 1WL(001) peak char- 
acteristic of the nano-layered structure of the crystallites is 
clearly visible. It is observed to be much more intense along 
the horizontal direction, which means that the directors of the 
crystallites are on average in a plane perpendicular to the direc- 
tion of the bundle (Fossum et ai, 2005). Furthermore, an axial 
symmetry around the horizontal direction perpendicular to the 
X-ray beam, that is, around the average direction for the bun- 
dles, is expected. Thus, the antinematic description presented 
in paragraph 2.2.2 is well suited to describe the orientational 
ordering in the chains. 

4.3.1. Method: The image analysis was done in the exact 
same way as for the nematic geometry studied previously. It is 
carried out on the ring corresponding to the first order reflection 
for the 1 WL hydration state, as shown by the arrow in Fig. 5(b). 
The choice of a Maier-Saupe functional form for the orientation 
distribution is not only phenomenological here. Indeed, this pre- 
cise shape is expected in a system of identical induced dipoles 
where alignment by an external field competes with the homo- 
geneizing entropy (Fossum et ai, 2005; Meheust et ai, 2005/?). 

4.3.2. Results: The normalized azimuthal profile extracted 
from the two-dimensional scattering image is shown in Fig. 10, 
with the corresponding model fitted on the whole range of 
azimuthal angles. Despite a significant noise away from the 
peak positions, the fit is satisfying. The uncertainties on the 
angular RMS width and order parameters (see table 2, column 
marked ER) are estimated in the way described in paragraph 
3.4.2. 

5. Discussion 
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The method reported in this paper displays several strong points 
with respect to a technique based on a one-dimensional record- 
ing, like the "rocking curve" technique. The first one is that it 
requires very little data acquisition time, since only one record- 
ing is necessary. Such a recording can be obtained in a couple 
of minutes at a 3rd generation synchrotron source. It therefore 
allows in situ measurements of the orientation distribution while 
changing experimental parameters. In the two examples given 
here, it could be the surrounding humidity and the magnitude 
of the applied electric field, respectively. The second benefit is 
that it allows determination of the reference orientation in the 
three dimensions of space. The rocking curve method provides 
the inclination ©o of that direction with respect to planes nor- 
mal to the incident beam in the case where 4>o = only. It 
does not intrinsically provide correct estimates of the angles ©o 
and 4>o without prior time-consuming adjustment of the sam- 
ple azimuthal orientation so as to ensure the condition <t>o = 0. 
Besides, it introduces a bias in the estimated width of the orien- 
tation distribution in the case where 4>o 7^ 0. 

We determine an uncertainty on the inferred ordering param- 
eters from the dispersion on the fitted values when changing 
the range of azimuth angles used for the fit. This uncertainty is 
therefore dependent on the noise level of the scattering image, 
and the process is more trustworthy if the diffraction peak 
considered is more intense. How the accuracy of the method 
depends on the peak's intensity with respect to the noise level 
is not quantified here. How the vicinity of another peak to the 
peak used to determine /, as in the OWL(OOl) peak in sec- 
tion 3, modifies the azimuthal profiles, is not known quanti- 
tatively either, and we propose our method as accurate when 
applied to a reflection whose peak is well separated from other 
peaks. Despite those shortcomings, we think the method has 
several advantages over the rocking curve method and allows to 
better constrain the description of the orientational order in the 
samples, because (i) the reference direction n is determined, (ii) 
no bias in the angular width is introduced in the case where n 
does not coincide with the Oz direction of the laboratory frame, 
(iii) an uncertainty is determined during the fitting procedure, 
for each azimuthal profile considered, and (iv) measurement 
consistency can be checked by comparing the results obtained 
from several different azimuthal profiles corresponding to sev- 
eral different diffraction peaks, on the same single recording. 
Finally, the method is easily carried out practically, since the 
only mathematical procedure involved is the fitting of Eq. 6 to 
azimuthal profiles. 

6. Conclusion and prospects 

We have developed a method to determine orientation distribu- 
tions in anisotropic powders of nano-layered crystallites from 
a single two-dimensional scattering image, in the case when 
it only depends on the deviation from a reference direction. 
The method relies on fitting the proper relation to the azimuthal 
decay of a given diffraction peak's amplitude. It allows determi- 
nation of the reference direction and of the angular dispersion of 
the distribution. It was applied succesfully to data obtained from 
two synchrotron X-ray scattering experiments, on two different 
systems of synthetic smectite clay particles where the particle 



assemblies exhibit a nematic and anti-nematic ordering, respec- 
tively. In the first case, the consistency of the inferred quantities 
was checked by comparing the results obtained from three dif- 
ferent orders of the same reflection. 

The method is promising as it is performed on a single two- 
dimensional recording, and therefore allows in situ determina- 
tion of orientational ordering in samples when thermodynami- 
cal conditions are varying. It shall be used in this respect, in the 
future. 



Appendix A 
Initial parameters for the fits 



The initial parameters for the fitting function in the form (15) 
are set to prior estimates of those parameters. The prior esti- 
mates are obtained from the data profiles using a few approxi- 
mations, as explained below. 

The parameter 1^ is first estimated as the minimum value 
of the data plot. A running average is performed on the data 
in order to smooth out the noise and obtain a curve with only 
two maxima: those corresponding to the peaks. The value for 
<£>q^ is set to the azimuthal position of the largest peak. Its peak 
intensity, I\ and that of the second peak, h are determined. The 
half-width at half maximum w of the largest peak is estimated 
numerically. The parameter is then defined from w by stat- 
ing: 

exp(mW) 



exp(W ! ' cos 2 w) 
which provides the relation 



m 



(0 



In 2 



1 — cos 2 w 



The maxima I\ and I2 are obtained for = <t>o and 1 
hence one can state: 



/i =C exp(m (,) cos 2 ' 



\0 - ®o')) 
h - = C exp (m» cos 2 (0 + 

From this, it follows that 



(16) 
(17) 

(18) 



In 



h 



= mW S in(20) sin(2©| ) '' ) ) 



1 n l 

from which an estimate of ©0 is obtained: 



(19) 



©W = - Asin 



In 



7 1 "^0 



mW sin(2fl) 

Finally, going back to (18), we obtain in the form 



(20) 



exp (mW cos 2 (19 - 0^)) 



(21) 
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Table 1 

Characteristic momentum transfer (in nm~ 1 ) and length scales (in A) measured 
for the OWL diffraction peaks from dry-pressed samples of Na-fluorohectorite. 
"Dry" denotes the dry-pressed Na-fluorohectorite samples, while "ER" denotes 
the electro-rheological bundles of Na-fluorohectorite particles. 



Reflection 


OWL(OOl) 


OWL(002) 


OWL(003) 


1 


6.50 ± 0.03 


13.02 ± 0.03 


19.57 ±0.03 


d 


9.67 ± 0.04 


9.65 ± 0.03 


9.63 ± 0.03 



Table 2 

Values obtained for the fit parameters m, ©o and 3>o, and for the angular width 
Wf and order parameter S. 



Profile 


Dry (001) 


Dry (002) 


Dry (003) 


ER (001) 


m 


10.45 ± 0.06 


10.21 ±0.27 


10.21 ±0.13 


3.10 ±0.23 


Wf 


19.10 ±0.07 


19.38 ±0.33 


19.37 ±0.16 


24.48 ± 0.50 


S 


0.85 ±0.01 


0.85 ±0.01 


0.84 ± 0.01 


0.56 ± 0.02 


00 


6.12 ±0.04 


5.57 ±0.10 


5.87 ± 0.05 


-1.64 ±0.05 


4>o 


176.37 ±0.01 


176.32 ±0.01 


176.37 ± 0.02 


267.74 ± 0.02 
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Figure 3 

Computed profiles describing how the intensity of a given diffraction peak (i,e, 
at a given deviation angle 2 9) depends on the azimuthal angle for a popula- 
tion of crystallites with a nematic orientational ordering. The computation was 
done for various values of the orientation (Qo.'&o) of the main director with 
respect to the fixed frame of the laboratory, and of the 2D RMS width, Wf, 
of the nematic angular distribution of the scatterers' directors. Two values of 
the deviation angle, 2 9 = 6.58°and 3.29°, were used, which means that two 
diffraction peaks corresponding to two different orders of the same reflection, 
were used. 



Figure 1 

In the fixed frame (x.y.z) attached to the laboratory (with x parallel to the inci- 
dent beam), the angular direction of the diffracted beam, denoted by point M in 
the detector plane yO'z, is defined by the deviation angle 2 9 and the azimuthal 
angle <f>. The area painted in light gray denotes the plane containing the incident 
and diffracted beams; it is normal to the yOz plane and makes an angle 4> with 
the z.Ox plane. The area painted in dark gray also denotes a plane normal to the 
to the yOz plane and that makes an angle with the zOx plane. It is used to 
define the orientation of the director n of a given crystallite with respect to the 
fixed frame, according to the angles © and <I>. 




Figure 2 

(a) The condition of specular reflection on the Bragg planes requests that 
= 9, as sketched here in the xOu plane, (b) The thick curve is that described 
on the unit sphere by the end of vector n (director of the scatterer) as the 
azimuthal position of the diffracted beam, <f>, varies between and 7r/2. Here the 
deviation angle 2 9 corresponds to the Bragg reflexion (002) on a nano-layered 
Na-Fluorohectorite particle in the hydration state 1WL, with A = 0.71 A. 




Figure 4 

Side view sketch of the experimental scattering geometry used to scatter X-ray 
by dry-pressed fluorohectorite samples. The insert shows a two-dimensional cut 
of the sample with the nano-layered crystallites lying on average aligned along 
a plane close to horizontal. 
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Figure 5 

Center part of two-dimensional scattering images obtained from the anisotropic 
powders of clay nano-stacks. The anisotropic rings used to infer the orientation 
distributions of the crystallites are denoted by arrows, (a) Picture obtained from 
dry-pressed samples. The camera length was 375 mm. (b) Picture from elec- 
trorheological chains of crystallites in oil. The camera length was 376 mm. The 
anisotropic ring in the center of (b) is the (001) ring for the nano-stack structure 
in the 1WL hydration state. Its anisotropy is orthogonal to that observed in (a). 
The image magnification is 1.5 with respect to (a). The broad peak outside the 
(001) ring is due to diffuse scattering from the silicon oil. 
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Figure 7 

Azimuthal profiles observed from the scattering of dry-pressed samples at three 
values of q corresponding to the three first orders of reflection by the population 
of crystallites in the 1WL hydration state. The profiles have been normalized 
and translated vertically for clarity. The fits according to Eq. (15) are shown in 
solid lines on the symbol plots for the data. 




Figure 8 

Close view of a chain formed under the application of a strong electric field 
by Na-fluorohectorite crystallites suspended in silicon oil. The original color 
image has been filtered. The gap between the electrodes, whose boundaries are 
shown by two vertical bold lines on the sides of the picture, is 2 mm. 




Figure 6 

Powder diffraction profiles recorded along different azimuthal directions on the 
2D detector, with a dry-pressed sample at T = 97.7°C. Several peaks character- 
istic of the nano-layered structure of the crystallites are visible. Their intensity 
decreases dramatatically as <j> increases due to the dropdown of the number of 
crystallites that meet the Bragg condition at these angles. 



Figure 9 

Top view sketch of the experimental scattering geometry used to scatter X-ray 
by electrorheological chains of Na-fluorohectorite crystallites. The electrodes 
are vertical, and the chains of aggregated crystallites along a direction which is 
close to perpendicular to the electrodes. The disorder in the chains is exagger- 
ated in this sketch. 
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Figure 10 

Azimuthal profile at q = 5.10 nm _1 , corresponding to the (001) peak of Na- 
fluorohectorite in the 1WL hydration state. The corresponding fit is plotted as a 
solid line on top of the data. 
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